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Search for pair production of first or second generation leptoquarks in proton— proton 
collisions at ^/s = 7 TeV using the ATLAS detector at the LHC 
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This paper describes searches for the pair production of first or second generation scalar lep- 
toquarks using 35 pb~^ of proton-proton collision data recorded by the ATLAS detector at 
y's = 7 TeV. Leptoquarks are seached in events with two oppositely charged muons or electrons 
and at least two jets, and in events with one muon or electron, missing transverse momentum and 
at least two jets. After event selection, the observed yields are consistent with the predicted back- 
grounds. Leptoquark production is excluded at the 95% CL for masses Mlq < 376 (319) GeV 
and Mlq < 422 (362) GeV for first and second generation scalar leptoquarks, respectively, when 
assuming the branching fraction of a leptoquark to a charged lepton is equal to 1.0 (0.5). 

PACS numbers: 14.80.Sv, 13.85.-t 



I. INTRODUCTION 

The Standard Model is extremely successful at describ- 
ing the elementary particles and their interactions, yet it 
gives no explanation for the striking symmetry between 
quarks {q) and leptons {£). This has, in part, been the 
motivation for many beyond-the-standard-model theo- 
ries that posit the existence of leptoquarks (LQ), parti- 
cles that carry both lepton and baryon quantum num- 
bers, which couple to both quarks and leptons, and 
carry the triplet charge of Quantum Chromo dynamics 
(QCD), color. Theories that predict leptoquarks include 
models that contain quark and lepton sub-structure 
theories that seek grand unification and models of 
extended technicolor j^]- Experimental bounds from 
searches for flavor-changing-neutral currents and lepton- 
family-number violation place restrictive limits on lepto- 
quark decays to different generations of quarks and lep- 
tons J3,^. Direct leptoquark searches at electron-proton 
colliders have sensitivity to first generation leptoquarks, 
but have a non-negligible dependence on the value of the 
LQ — £ — q coupling [(^ , whereas second and third genera- 
tion leptoquarks are produced via lepton-flavor-violation 
mechanisms Q- The relatively large cross sections [S] 
for scalar leptoquark production from proton-proton col- 
lisions lead to the expectation that early LHC data offer 
sensitivity to a mass range beyond that probed by other 
accelerators. Currently the most stringent limits come 
from the Tevatron [l^ and from the CMS experiment 
at the LHC [Hllil. 

Scalar leptoquarks can be produced in proton- 
proton collisions either in leptoquark-antileptoquark 
pairs (LQ LQ) or singly. Single LQ production involves 
an unknown LQ — i — q coupling, whereas the pair pro- 
duction reaction occurs mostly via QCD processes which 
involve only the strong coupling constant. Therefore, 
results from LQLQ searches have negligible model de- 
pendence, and the only relevant parameter for scalar 
LQ production is the leptoquark mass Mlq [13|. Final 



states from LQ pair production have either two same- 
flavor oppositely-charged leptons and two jets {Ujj)', a 
lepton, a neutrino and two jets (Ivjj)] or two neutrinos 
and two jets. The leptons / here and throughout the pa- 
per are either electrons for first generation LQ or muons 
for the second generation [ij . The cross sections for the 
leptoquark-mediated processes pp — )• Ujj and Ivjj can 
be written as ctlq x Z?^ and ctlq x 2/3(1 — respectively, 
where /? is the branching fraction for a single leptoquark 
to decay into a charged lepton and a quark. 



II. ANALYSIS STRATEGY 

This paper reports a search for scalar leptoquark pair 
production carried out using 35 pb^-'^ of data recorded by 
the ATLAS detector during the 2010 LHC proton-proton 
running period. The analysis is performed separately in 
the Ujj final state and in the liyjj final state. These 
searches are combined, leading to final results presented 
as a function of j3 and M^q for the first and second gen- 
erations. 

Analyses for both final states begin by selecting event 
samples that have high acceptance for signal production. 
At this initial stage, these samples are dominated by the 
major backgrounds, Z-|-jets and ti for the Ujj case, and 
M^-|-jets and tt for the Ivjj case. The samples are then 
subdivided into orthogonal control and signal regions. 

The control regions are used to validate the back- 
ground modeling by the Monte Carlo (MC) simulation, 
and the signal region is used to search for evidence of LQ 
production. The signal region is defined using an a priori 
optimization procedure based on simulated background 
and signal events. 

A major difference between LQ events and back- 
grounds is the presence of jet-lepton (also jet-neutrino 
in the case of Ivjj) pairs coming from the decay of 
the parent LQ, giving a peak in the reconstructed jet- 
lepton mass spectrum for the signal. The approximate 
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reconstruction of these masses provides the most impor- 
tant variables used to distinguish signal and background 
events. In addition, large LQ masses give rise to larger to- 
tal measured transverse energy in LQ pair events than is 
seen for background events, giving another means to dis- 
tinguish signal from background. Finally, reconstructed 
boson masses can be used to reject the dominant back- 
grounds from V{V = W^, 2')-|-jets production where the 
boson decays into leptons. 



In the lljj channel, an average reconstructed lepto- 
quark mass M lq is defined for each event by com- 
puting the average of the masses from the two lepton- 
jet combinations in the event. Both possible assign- 
ments of the two leptons and the two leading jets to 
LQ parents are considered. The chosen assignment is 
that which gives the smallest absolute difference between 
the two reconstructed masses. The probability to get 
the correct pairing with this method is of the order of 
90%. The transverse energy in an event is defined as 
the scalar sum of the transverse energy (momentum) of 
the two electrons (muons) and of the two leading jets, 
— + pf^ + pl^ + pl^ , where the transverse plane is 
defined as relative to the beam axis The invariant 
mass of the dilepton pair AIu provides rejection of the 
Z+^ets background. 



III. ATLAS DETECTOR 

ATLAS [l^ [13 is a multi-purpose detector with a 
forward-backward symmetric cylindrical geometry and 
nearly An coverage in solid angle. The three major 
sub-components of ATLAS are the tracking detector, 
the calorimeter, and the muon spectrometer. Charged 
particle tracks and vertices are reconstructed in silicon 
based tracking detectors that cover |r/| < 2.5 and tran- 
sition radiation detectors extending to I77I < 2.0. The 
inner tracking system is immersed in a homogeneous 2 T 
magnetic field provided by a solenoid. Electron, photon 
and jet transverse energies are measured in the calorime- 
ter. The ATLAS calorimeter system is segmented into 
a central barrel and endcaps collectively covering the 
pseudorapidity range I77I < 4.9, and is equipped with 
both electromagnetic and hadronic calorimeters. Sur- 
rounding the calorimeter, a muon spectrometer with air 
core toroids, a system of precision tracking chambers, 
and detectors with triggering capabilities provides precise 
muon identification and measurements. A three-level 
event-triggering system allows for the selection of inclu- 
sive muon events with muon transverse momentum (pl^) 
greater than 13 GeV and inclusive electron events with 
electron transverse energy (Et) greater than 15 GeV. 



IV. SIMULATED SAMPLES 



In the Ivjj channel, LQ mass equivalents are also 
defined. The neutrino transverse momentum p!^ is in- 
ferred from the missing transverse momentum in the 
event E^^^^ as described in Section V. As in the lljj 
final state, two pairings of lepton and jet and E^^^^ 
and jet are possible. However, because the component 
of the neutrino momentum along the beamline is un- 
determined, only one mass Mlq from the charged lep- 
ton and a jet can be reconstructed. The E^'^^ and the 
remaining jet are used to compute a transverse mass, 

MjJq = \J 2p^rj.E^^^^ {1 — cos 4^ ) in which pl^ is the trans- 
verse momentum of a jet, and is the angle between 
the plj, and the E™^^^ vectors in the transverse plane. 
In analogy with the lljj final state, the chosen pairing 
is that which gives the smallest absolute difference be- 
tween A/lq and M^q, also resulting in the correct as- 
signment more than 90% of the time. The measured 
transverse energy is defined as the scalar sum of the lep- 
ton transverse momentum, the missing transverse mo- 
mentum and the momentum of the two leading jets in 
the event, 5^ = pf^ + -B^"'''' +j^r^+j^^. An additional 

transverse mass variable Mt — \J 2p^^i?™'*'^(l — cos (j)^) 
provides rejection against the dominant W^-l-jets back- 
ground. Here, p^ is the measured lepton transverse mo- 
mentum, and (f/ is the angle between the p^ and E"™'^** 
vectors. 



Event samples processed through a detailed detector 
simulation [18,] based on GEANT4 [19i] are used to de- 
termine most background and all signal yields. As the 
recorded data contain a non-negligible fraction of events 
with more than one interaction per beam crossing, sim- 
ulated events passing the initial lepton selection criteria 
are reweighted in order to bring the vertex multiplicity 
distributions in data and simulation into agreement. All 
background samples use ATLAS tune MC09 [1^. The 
vector boson processes, Vl^-l-jets and Z-|-jets, are gen- 
erated using ALPGEN v2.13 [IH with CTEQ6L1 [H 
for the parton distribution functions (PDFs), interfaced 
to HERWIG V6.510 23] and JIMMY v4.31 ]24]- Ex- 
clusive samples with zero to four additional partons (np 
with n = 0, 1, 2, 3, 4) and an inclusive sample with five 
or more additional partons (np with n = 5) are used. 
The cross sections are computed using the ALPGEN 
cross sections scaled so that the sum of the rtp sample 
cross sections are equal to the next-to-next-to-leading- 
order (NNLO) inclusive cross sections times branching 
fraction to a single lepton species: a{W — > Iv) = 
10.46 ± 0.42 nb and (jI^/tI ^ U) = 1.070 ± 0.054 nb 
for Mil > 40 GeV [21, Additional samples for 

cross checks are generated using SHERPA 1.1.3 27] with 
PDF set CTEQ6L1, and PYTHIA 6.421 ^ with PDF 
MRST 2007 LO* ^29,]. Samples of tt events are produced 
using MC@NLO v3.41 [30| and POWHEG 1.01 patch 
4 ^ with PDF CTEQ6.6M [H interfaced to HER- 
WIG for parton showering. Cross-sections are scaled 
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to the next-to-next-to-leading-log (NNLL) prediction of 
165l^g pb m, [S^. Single top events are generated us- 
ing MC@NLO with cross sections of 3.94 pb, 58.7 pb 
and 13.1 pb for the s-, t- and tW^-channcls respec- 
tively. Their uncertainties are ±10% [35„^]. Diboson 
events are generated using HERWIG. Next-to-leadin g or - 
der (NLO) cross sections are calculated with MCFM [37i |: 
44.9 ±2.2 pb, 18.0 ±1.3 pb, and 5.96 ±0.30 pb for WW, 
WZ (M« > 40 GeV) and ZZ {Mu > 60 GeV) respec- 
tively. Cross section uncertainties take into account scale 
and PDF uncertainties, as well as differences with other 
generators. 

Signal events for LQ masses of 250 to 400 GeV with 
a 50 GeV binning are generated with PYTHIA and tune 
D6 [38] with cross sections and uncertainties determined 
from Ref. ^S] using the CTEQ6.6 ^39,] PDF set. The un- 
known LQ — f—g coupling value is set to 0.01 x y/ Ana em- 
This corresponds to an LQ full width of less than 1 MeV, 
and a negligible decay length i40, 41j. 



V. OBJECT IDENTIFICATION 

Electron candidates are reconstructed from energy 
clusters in the electromagnetic calorimeter. Electron 
identification 42] is performed using the transverse 
shower shape, the longitudinal leakage into the hadronic 
calorimeter, and the requirement that a good-quality 
track points to the cluster. The electron transverse en- 
ergy is measured in the calorimeter, while its direction 
is obtained from the track. Further rejection against 
hadrons is achieved by using the energy deposit patterns 
in the first layer of the EM calorimeter. In order to sup- 
press the background from photon conversions, a hit in 
the first layer of the pixel detector is required. Electrons 
used in this analysis are required to have Et > 20 GeV 
and \ri\ < 2.47, with the exclusion of the poorly in- 
strumented region between the barrel and the end-cap 
calorimeters at 1.35< jry] <1.52. A small fraction of 
events with electrons near problematic regions of the elec- 
tromagnetic calorimeter readout is removed. 

Muons selected for this analysis are required to have 
pt^ > 20 GeV and \ri\ < 2.4. Muon tracks are recon- 
structed independently in the inner detector and in the 
muon spectrometer, with a minimum number of hits re- 
quired in each. A good match is required between the 
tracks found in the inner tracker and the muon spectro- 
meter. In order to reject the cosmic ray background, 
tracks from muon candidates must extrapolate back to 
the reconstructed event vertex, satisfying \do\ < 0.1 mm 
and \za\ < 1 cm, where do is the minimum distance be- 
tween the muon trajectory and the event primary vertex 
in the plane perpendicular to the beam direction, and zq 
is the corresponding distance parallel to the beam direc- 
tion. 

Finally, both electrons and muons are required to be 
isolated from other energy in the calorimeters by im- 
posing E^°'"'/E'T < 0.2 and Ef'^^/pt^ < 0.25, for elec- 



trons and muons, respectively. Here E'™"*^ is the trans- 
verse energy in the calorimeter in a cone of size AR = 
y/ {Ai]y + (A0)2 — 0.2 centered on the lepton direction, 
excluding the lepton contribution. 

The presence of neutrinos is inferred from the missing 
transverse momentum E^^^^. The E!^^^^ is defined from 
the vector sum of the transverse energy in calorimeter 
cells included in topological clusters [la] and the trans- 

cells 

verse momentum of the muon, E^^^^ = — (ESt +Pt)- 
The clusters are corrected to take into account the dif- 
ferent response to hadrons compared to electrons or pho- 
tons, as well as dead material and out-of-cluster energy 
losses HI]. 

Jets are reconstructed from calorimeter energy clusters 
using the anti-fct [3 ] algorithm with a radius parameter 
R — 0.4. After applying quality requirements based on 
shower shape and signal timing with respect to the beam 
crossing [4^, jets selected for this analysis must satisfy 
PT > 20 GeV and \ri\ < 2.8 and must be separated from 
leptons by AR > 0.5. Galibrations of lepton and jet 
transverse momenta, which are mostly derived from con- 
trol samples in data, are applied prior to making the 
kinematic selections. 



VI. PRESELECTION 

Initial selection criteria define event samples with high 
signal acceptance and yields which are dominated by 
the major backgrounds F+jets and tt. Events in both 
electron and muon analyses are required to have at 
least one reconstructed proton-proton interaction ver- 
tex with at least three charged-particle tracks associ- 
ated to it. The dilepton channels additionally require 
exactly two electrons (muons) with £'t(p!^) > 20 GeV. 
The single lepton channels require exactly one electron 
(muon) with Eripi^^) > 20 GeV, i;^'"*^ > 25 GeV and 
Mt > 40 GeV. For the evjj channel only, a triangle cut 
A0(jet,£'^'"") >1.5x(l-i;™''745), where </> is in radians 
and E^^^^ in GeV is applied, in order to reduce the mul- 
tijet background contamination. All analyses require at 
least two jets satisfying pT > 20 GeV. The two highest 
Pt jets are used to calculate LQ masses. 

The acceptance for both signal and background is es- 
timated by applying all selection criteria to simulated 
events. Expected yields are obtained by scaling the ac- 
ceptance by the predicted cross section and integrated 
luminosity. 

After preselection, good agreement is observed be- 
tween simulation and data in both the number of events 
and the shape of distributions of kinematic variables. 
The determination of the expected backgrounds is dis- 
cussed in Sections VII and VIII. The expected and ob- 
served yield for the four different channels are summa- 
rized in Table m 
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TABLE I: The predicted and observed yields for the prese- 
lected sample for all channels. Both statistical and systematic 
uncertainties are included. 



Channel 


Predicted Yield 


Observed Yield 




610±240 
610011?* 
8301?* 
9500±2500 


626 
6088 
853 
9248 



VII. BACKGROUND DETERMINATION 

Small differences between data and simulation for reso- 
lutions and trigger and reconstruction efficiencies are de- 
termined in control data samples and applied to the sim- 
ulated events. These corrections influence the obtained 
yields by less than 2%. The background arising from 
lepton misidentification is determined using data-driven 
methods in the four channels, as is the contribution from 
the Z+jets background to lljj final states. All other 
backgrounds are modeled with Monte Carlo simulations 
and tested with data in control regions defined to enhance 
their contributions, as described in Section VIII. 

Because the Z+}ets yield in the signal region for the 
lljj final state arises in the tails of distributions, it is 
determined using a data-assisted method. The yield in 
the signal region A^^* is calculated as: 

MC 

where A^ and -/Vf/^ are the observed numbers of events 
in data and MC, respectively, in a 20 GeV wide dilep- 
ton mass window around the nominal Z boson mass, 
and N^^jfj is the expected Z-f jets contribution to the sig- 
nal region (defined in Table irv|) . To estimate a system- 
atic uncertainty, the prediction of N^j^ was derived with 
and without requirements on the number of jets in the 
event, as well as with the different Monte Carlo genera- 
tors described in Section IV. The 3% background in the 
20 GeV mass window from ti and diboson events, esti- 
mated with Monte Carlo, is subtracted from the data, 
as is the <1% contamination from fake leptons. The 
data-assisted methods (with and without a cut on the 
number of jets and with various mass window definitions) 
and the purely Monte Carlo-based estimates agree within 
10%. The largest difference is observed between results 
obtained with ALPGEN and SHERPA, and is used as 
a systematic uncertainty. The final estimate is obtained 
using ALPGEN for Nfjc and A^'/^. 

Jets misidentified as leptons and leptons from semilep- 
tonic decays of hadrons are referred to as fake leptons. 
In this paper, events with fake leptons are referred to 
as QCD background [isj. The QCD background in the 
lljj channel is estimated using a fitting method. Tem- 
plates for both real and fake electrons (muons) are made 
for the E^'^^/Et (£'™"°/pt) variable, before applying 



the isolation cut. Simulation is used to determine tem- 
plates for real leptons. Templates for fake leptons are 
derived from the data by selecting events with exactly 
one lepton. Events containing W candidates are rejected 
by requiring E^^^^ < 10 GeV. The small residual real 
lepton contamination is estimated from simulated events 
and is subtracted from the data. In the eejj channel, 
additional templates are made for the fraction of hits 
associated with transition radiation from electrons. For 
both channels, fits with these templates are performed 
on dilepton events with two or more jets giving the prob- 
ability that an event contains at least one fake lepton. 
Fits are made inside the Z mass window, as well as in 
the signal region. Both fits indicate a small expected 
contribution from fake leptons. In the signal region, an 
independent upper limit on the number of fakes is set 
by extrapolating from the number of fakes determined in 
lower jet multiplicity bins and with lower dilepton mass. 
The QCD background contribution to the fJ-vjj analysis 
is determined using a scaled control sample method. In 
this method, a pair of variables that cleanly separate real 
muons from vector boson or leptoquark decay from fake 
muons is chosen to divide the full sample into four statis- 
tically independent regions. One of the four regions cor- 
responds to the signal region, which also contains muons 
from vector-boson decays. The background in the signal 
region is then found by scaling the yield in one of the 
non-signal regions by the ratio of yields in the remaining 
two regions. This analysis uses E^^^^ and muon impact 
parameter \do\ to define the different regions. The three 
background regions are defined as: |cZo| < 0.1 mm and 
E'j}''''' < 25 GeV (A), \do\ > 0.1 mm and E!^''''' < 25 GeV 
(B) and Idol > 0.1 mm and > 25 GeV (C). The 

signal region (D) is the region with |do| < 0.1 mm and 
^miss ^ 25 GeV. The correlation between these two vari- 
ables is less than 10% in simulated samples, giving sta- 
bility to the method. The yield Ad in the signal region 
is given by Ad = ^j^;^* , where A refer to the observed 
numbers of events in the four statistically independent 
regions after subtraction of the small, residual contri- 
butions from events with real muons in the background 
regions, which are estimated using simulated events. 

In the evjj channel, the normalization of the QCD 
background is determined from a fit to data of the sum 
of the Mr distributions for QCD events and all other 
backgrounds, primarily M^-f jets and tt. This sum is con- 
strained to equal the total data yield, with the QCD frac- 
tion being the fit parameter. The template for the QCD 
background shape is determined using a QCD enriched 
sample constructed by taking the difference in shapes 
between two samples, a loose sample selected using only 
the electrons passing the trigger requirements, and the 
nominal sample, selected using the full electron identifi- 
cation requirements. Residual contamination from real 
electrons in the QCD enriched sample is estimated to be 
7%. The latter is estimated with a loose-tight matrix 
method [s^, and used to perform a shape-dependent 
subtraction. For both muon channels, the background 
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from events containing cosmic ray muons is negligible. 



VIII. CONTROL REGIONS 

Control samples are used to validate the background 
determination based on MC. The two most important 
control samples for the lljj analyses are (1) Z+jet: 
events for which the dilepton invariant mass lies within 
the Z mass window 81 < Ma < 101 GeV and (2) tt: 
events that contain at least two jets and both an electron 
and a muon selected, as defined in Section V. By defini- 
tion, the tt control region is common for both the eejj 
and fifijj channels. The most important control samples 
for the Ivjj analyses are (1) W+2 jets: events with ex- 
actly two jets, a charged lepton and E^^^^ such that Mt 
is in the region of the W Jacobian peak (40 GeV < A/t < 
150 GeV), (2) W+3 jets: as in (1) but with at least 3 jets 
and (3) a ti enriched sample which requires at least four 
jets with > 50 GeV, pi^ > 40 GeV and > 30 GeV. 
The expected signal contamination in the control regions 
is at most 1%. The predicted and observed yields in these 
control samples are shown in Tables and IIIIl Distri- 
butions of are shown in Fig. [1] for the lljj samples 
in the Z+^et and ti control regions, and distributions for 
the reconstructed Mlq are shown in Fig. [5] for the li'jj 
samples in the 1^-1-2 jets and the ti control regions. 



IX. SELECTION OPTIMIZATION 

The signal regions in the four final states are defined 
using a random grid search optimization method [47j . 
The optimization is performed using simulated signal 
and background samples (determined as described in Sec- 
tions VII and VIII) after preselection, with the final cuts 
determined by minimizing the Poisson probability that 
the predicted background fluctuates to at least the pre- 
dicted signal plus background. A 10% overall scale un- 
certainty on the yields, roughly matching the uncertainty 
on the luminosity, is included in the optimization as a 
systematic uncertainty. Optimization for the single lep- 
ton (dilepton) channels nominally is done using signal 
samples with LQ masses of 300 (400) GeV. Cuts derived 
using different LQ masses typically differ from the nom- 
inal values within statistical uncertainties. A variety of 
input variables are considered. The requirements that 
give the best separation between signal and background 
are shown in Table IIVI The yields obtained after final 
selection for the various channels are shown in Table fVl 
along with the predicted background and signal. 



X. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are derived for a variety of 
sources, including data-simulation differences in trigger 



and reconstruction efficiencies and in the energy and mo- 
mentum resolutions for leptons, jets and E^^^^, instanta- 
neous and integrated luminosity, modeling of the under- 
lying event, variations in the method used to determine 
the QCD background, and uncertainty on the LQ pair- 
production cross section. A summary of the systematic 
uncertainties for the four final states is shown in Table IVTl 

The lepton trigger and reconstruction efficiency sys- 
tematic uncertainties are derived by varying the selection 
of the Z event sample used to measure in situ efficien- 
cies, and by varying the treatment of the (small) back- 
ground in this sample. In addition, a 1% uncertainty is 
included on the muon isolation requirement that accounts 
for the difference of the efficiency of the isolation require- 
ment in data and ti and LQ simulated samples. Finally, 
a 3% uncertainty is added to account for differences in 
the muon \do\ distributions between data and simulation. 
Lepton momentum scale and resolution uncertainties are 
obtained by comparing the peak and width of Z ^ i£ 
events in data and Monte Carlo. 

The jet energy scale and resolution are varied by their 
uncertainties [H, |49i] for all simulated events, and their 
impact estimated independently. In addition a 5% un- 
certainty is added in quadrature to the jet energy scale 
uncertainties to account for differences in response for 
quark and gluon jets. These variations in scale and res- 
olution are also propagated to the E^^^^^. The system- 
atic uncertainty from instantaneous luminosity effects is 
evaluated by comparing the results from simulated sam- 
ples with and without additional minimum bias events 
(pileup) added. This is an overestimate of the system- 
atic, but is still small (2-6%, depending on the sample) 
compared to other sources. 

Systematic uncertainties on the QCD background are 
determined by comparing results from alternate normal- 
izations to those described in Section VII. The uncer- 
tainty on the lljj final states is determined by varying 
the estimate between the nominal prediction and the up- 
per level from the extrapolation method described in Sec- 
tion VII. The uncertainty in the evjj final state is deter- 
mined by comparing the default method, based on fits to 
the Mt distribution, to alternate fits using the E!^^^^ or 
the electron Et distributions. The largest fractional dif- 
ference (22%) between the nominal fit and fits using the 
alternate variables is taken as the systematic uncertainty. 
The uncertainty in the background to the fii'jj final state 
is determined by comparing the default method to one 
that uses the muon isolation instead of the c?o variable 
The difference in yield between the two estimates is used 
to determine the systematic uncertainty, giving 27%. 

The systematic uncertainties for the production mod- 
els of iy-|-jets and Z-|-jets events in the single lepton 
analysis are determined by varying parameters at the 
generator-level as described in [26|, including variations 
of the renormalization and factorization scales and min- 
imum parton px- The largest absolute deviation (40%) 
is used as a systematic uncertainty. In the dilepton anal- 
ysis, the systematic uncertainty is determined by vary- 
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TABLE II: The predicted and observed yields in the control samples for the electron final states. Top refers to both single top 
and ti events. Both statistical and systematic uncertainties are included. 





eejj 




Event 


Control Region 




Control Region 






Source 


Z+>2 jets 




ti 




W + 2 jets 


W+ > 3 jet 




tt 




y+jets 


150 ± 23 


0.3 


± 


0.1 


2100 ± 700 


580 ± 190 


180 


± 


60 


Top 


2.0 ± 0.3 


24 


± 


4 


21 ± 4 


44 ± 9 


210 


± 


40 


Diboson 


2.0 ± 0.3 


0.8 


± 


0.1 


17 ± 4 


8.3 ± 1.9 


2.1 


± 


0.5 


QCD 


4.0 t i*o° 


0.0 


+ 


0.1 
0.0 


64 ± 14 


68 ± 15 


29 


± 


7 


Total Bkg 


158 ± 25 


25 


± 


4 


2200 ± 700 


700 ± 200 


420 


± 


80 


Data 


140 




22 




2344 


722 


425 





TABLE III: The predicted and observed yields in the control samples for the niuon final states. Top refers to both single top 
and ti events. Both statistical and systematic uncertainties are included. 









Event 


Control Region 




Control Region 






Source 


Z+>2 jets 




ti 




W + 2 jets 


W+ > 3 jet 


ti 




V+jets 


190 ± 24 


0.3 


± 


0.1 


3300 ± 1100 


900 ± 300 


250 ± 


80 


Top 


2.7 ± 0.5 


24 


± 


4 


14 ± 3 


53 ± 1 


260 ± 


50 


Diboson 


0.2 ± 0.1 


0.8 


± 


0.1 


28 ± 6 


14 ± 3 


3.0 ± 


0.7 


QCD 


6.0 t 


0.0 


+ 


0.1 
0.0 


300 ± 100 


130 ± 50 


54 ± 


32 


Total Bkg 


200 ± 25 


25 


± 


4 


3600 ± 1100 


1100 ± 330 


570 ± 


120 


Data 


216 




22 




3588 


1120 


547 



ing the dilepton mass window and replacing the ALP- 
GEN event generator with SHERPA. The difference be- 
tween ALPGEN and SHERPA dominates, giving an un- 
certainty in the eejj{fifijj) channel of 34 (45)% on the 
Z+jets prediction in the signal region. 

Systematic uncertainties are evaluated for the ti pro- 
duction model by comparing the predictions obtained 
with MC@NLO and POWHEG. The resuh is a 35% un- 
certainty for both Ujj and liyjj channels. 

An integrated luminosity uncertainty of 11% ^50] is 
applied to all backgrounds determined from simulated 
events and to the signal. Additional signal uncertainties 
are obtained from varying the renormalization scale pa- 
rameter to 0.5AfLQ and to 2AfLQ (15%), from different 
PDF choices (13-17%) for LQ masses in the range 300- 
400 GeV, and from Initial- and Final-State- Radiation ef- 
fects (2%). The systematic uncertainties on the theory 
cross sections used to normalize backgrounds are given 
in Section IV. 



XI. RESULTS 

The expected and observed numbers of events are 
shown in Table |Vl Data and Standard Model expec- 
tations are in good agreement. 95% CL upper limits on 
LQ pair-production cross sections are determined using 
a modified frequentist approach [sil [52j . Systematic un- 
certainties are incorporated in the limit calculation as 
nuisance parameters and are integrated out. The lim- 



its are translated into bounds in the /? versus LQ mass 
plane. 

Monte Carlo studies show that a better sensitivity is 
achieved when using kinematic shapes rather than just 
the total number of events. For the eejj and fJ-fJ-jj fi- 
nal states, the observed and predicted St distributions 
are used in the limit-setting procedure, and for the eiyjj 
and ^vjj final states the observed and predicted Mlq 
distributions are used. The St distributions are shown 
in Fig. [3] for the Ujj channels. Figure 2] shows the Mlq 
distributions for the single lepton final states. 

The 95% GL upper bounds on the cross section for first 
and second generation LQ pair production as a function 
of mass are shown in Fig. [5] for /3 —1.0 and f3 —0.5. The 
combined limits are also shown in the /3 vs. Mlq plane in 
Fig. [HI for both generations. The expected and observed 
combined limits including all systematic uncertainties are 
shown in Table EII] for both (3 =1.0 and 0.5. The sys- 
tematic uncertainties, of the order of 50%, only change 
the limits by 5 to 10 GeV, depending on the value of (3. 



XII. CONCLUSIONS 

This paper reports the results of searches for pair pro- 
duction of first or second generation scalar leptoquarks 
using a data sample corresponding to an integrated lu- 
minosity of 35 pb""'^. The data in the high signal-to- 
background signal region are in good agreement with the 
Standard Model expectations. 95% GL upper bounds 
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on the production cross section are determined. These 
are translated into lower bounds for first (second) gener- 
ation leptoquark masses of Mlq > 376 (422) GeV and 
Mlq > 319 (362) GeV for /3 = 1.0 and (3 = 0.5, re- 
spectively. These are the most stringent bounds to date 
from direct searches for leptoquarks in much of the phase 
space. 
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FIG. 1: St distributions for the eejj (top) and the fJ.fJ.jj 
final states (middle) in the Z control region. The bottom 
figure shows the distribution for the dilepton tt control re- 
gion. The data are indicated by the points and the Standard 
Model backgrounds are shown with cumulative distributions. 
The QCD background is estimated from data, while the other 
background contributions are obtained from simulated sam- 
ples. The top background includes both tt and single top 
events. The expected contribution from a potential LQ sig- 
nal is negligible and is not included in the figure. The MC 
uncertainty band shows the combined statistical and system- 
atic uncertainties on the prediction in each bin. 
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FIG. 2: Reconstructed Mlq distributions for the ei/jj (top) and the fivjj (bottom) final states in the W + 2 jets (left) and tt 
(right) control regions. The data are indicated by the points and the Standard Model backgrounds are shown with cumulative 
distributions. The QCD background is estimated from data, while the other background contributions are obtained from 
simulated samples. The top background includes both ti and single top events. The expected contribution from a potential 
LQ signal is negligible and is not included in the figure. The MC uncertainty band indicates the combined statistical and 
systematic uncertainties on the prediction in each bin. 



TABLE IV: The selection requirements used to define the signal region, as obtained from the optimization procedure, > 
30 GeV implies that the pt of both the two leading leptons and the two leading jets in the dilepton samples should exceed 30 
GeV. 



eejj and fifijj 


evjj 




Mil > 120 GeV 
Mlq > 150 GeV 
pi|J' > 30 GeV 
S^T^ > 450 GeV 


Mt > 200 GeV 
Mlq > 180 GeV 
M^Q > 180 GeV 

S!^ > 410 GeV 


Mt > 160 GeV 
Mlq > 150 GeV 
M^fq > 150 GeV 

S!i > 400 GeV 
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TABLE V; The predicted and observed yields in the signal region for all channels. The lljj (li^jj) channel signal yields are 
computed assuming /3 = 1.0 (0.5). Both statistical and systematic uncertainties are included. Top refers to both single top and 
tt events. 



Source 


eejj evjj n^jj ^.vjj 


l/+jets 
Top 

Diboson 
QCD 


0.50 ± 0.28 
0.51 ± 0.23 
0.03 ± 0.01 
0.02 t 


0.65 ± 0.38 
0.67 ± 0.39 
0.10 ± 0.03 
0.06 ± 0.01 


0.28 ± 0.22 
0.52 ± 0.23 
0.04 ± 0.01 

0.00 1 g:SJJ 


2.6 ± 1.4 
1.6 ± 0.9 
0.10 ± 0.03 
0.0 ± 0.0 


Total Bkg 

Data 


1.1 ± 0.4 
2 


1.4 ± 0.5 
2 


0.8 ± 0.3 



4.4 ± 1.9 

4 


LQ(250 GeV) 
LQ(300 GeV) 
LQ(350 GeV) 
LQ(400 GeV) 


38 ± 8 
17 ± 4 
7.7 ± 1.7 
3.5 ± 0.8 


9.6 ± 2.1 
5.1 ± 1.1 
2.6 ± 0.6 


45 ± 10 
21 ± 5 
9.4 ± 2.1 
4.4 ± 1.0 


13 ± 3 
6.4 ± 1.4 
3.0 ± 0.7 



TABLE VL Systematic uncertainties for the lljj and the lujj final states. The lepton trigger, identification and momentum 
(energy) scale and resolution uncertainties are small and grouped together. Single top and events are also grouped together. 
Uncertainties marked with * (**) are only for the electron (muon) channels. QCD modeling systematics are not shown. All 
numbers are in percentages. 



Channel 


y+jets 

^iji ifji 


Top 
lljj 


>-^3j 


Diboson 


LQ (300 GoV) 
lljj It^jj 


Production Cross Section 




4 


13 


13 


5 


5 


18 


18 


Modeling 


34*, 45** 


40 


35 


35 










Electron Energy Scale & Resolution* 


+13, -0.2 


5 


10 


2 


7 


1 


8 


1 


Muon Momentum Scale &; Resolution** 


20 


5 


7 


2 


8 


1 


6.7 


1 


Jet Energy Scale 


6 


+22, -13 


+9, -18 


32 


+16, -6 


+ 17, -24 


2 


3 


Jet Energy Resolution 


16 


10 


0.3 


26 


4 


14 


0.3 


3 


Luminosity 


0.3 


11 


11 


11 


11 


11 


11 


11 


Pile up 


< 0.1 


5 


< 0.1 


4 


< 0.1 


6 


< 0.1 


2 


Total Systematics 


39* 
52** 


+49, -45 


47* 
(+49 -44)** 


57 


(+22, -16) +26,- 31 


22 


22 
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FIG. 3: St distribution for the eejj (left) and the /i/ijj final states (right) after all selections. The data are indicated by the 
points and the Standard Model backgrounds are shown with cumulative distributions. The expected LQ signals for various 
masses are also shown. 
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FIG. 4: Mlq distribution for the ei/jj (left) and the jJ.vjj final state (right) after all selections. The data are indicated by the 
points and the Standard Model backgrounds are shown with cumulative distributions. The expected LQ signals for various 
masses are also shown. 



TABLE VIL 95% CL expected and observed limits on the first and second generation LQs under different assumptions on /3. 
The expected limits include systematic uncertainties. 



Type (13) 


Expected limit (GeV) 


Observed limit (GeV) 


1st 


generation (1.0) 


387 


376 


1st 


generation (0.5) 


348 


319 


2nd 


generation (1.0) 


393 


422 


2nd 


generation (0.5) 


353 


362 
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FIG. 5: 95% CL upper bounds on the LQ pair production cross sections as a function of LQ mass for the combination of the 
eejj and the eujj channels (top), and for the fifijj and the fivjj channels (bottom). The figures on the left (right) show the 
limit for /3 =1.0 (0.5). The expected limit is indicated by the dashed line. The light yellow (dark green) solid band contains 
68% (95%) of possible outcomes from pseudo-experiments in which the yield is Poisson-fluctuated around the background-only 
expectation. Systematics uncertainties are included. The observed limit is indicated by the solid line. The theory prediction is 
indicated by the dotted line, which includes systematic uncertainties due to the choice of the PDF and to the renormalization 
and factorization scales (A(t,j). 
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FIG. 6: 95% CL exclusion region obtained from the combination of the two electron channels (left) and the muon channels 
(right) shown in the j3 versus leptoquark mass plane. The gray area indicates the DO exclusion limit [3, [13], and the thick 
dotted line the CMS exclusion [HI. The dotted and dotted-dashed lines show the individual limits for the II jj and the 
Ivjj channels, respectively. The combined expected limit is indicated by the thick dashed line. The solid band contains 68% of 
possible outcomes from pseudo-experiments in which the yield is Poisson-fluctuated around the background-only expectation. 
Systematic uncertainties are included. The combined observed limit is indicated by the solid line. 
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